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Abstract 

We study the Dirac equation with Coulomb-type vector and scalar 
potentials in D + 1 dimensions from an su(l, 1) algebraic approach. 
The generators of this algebra are constructed by using the Schrodinger 
factorization. The theory of unitary representations for the su(l, 1) 
Lie algebra allows us to obtain the energy spectrum and the supersym- 
metric ground state. For the cases where there exists either scalar or 
vector potential our results are reduced to those obtained by analytical 
techniques. 
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1 Introduction 

As it is well known the factorization methods are very important to study 
quantum systems [U [21 E] , since they are the basis for obtaining the energy 
spectrum and the eigenfunctions in an algebraic way [3]. 

The fundamental ideas of factorization in quantum physics were settled 
by Dirac [JJ and Schrodinger [2J. However, Infeld and Hull [3] were the first 
to introduce a systematic method to factorize and classify a large class of 
potentials. Moreover, for several problems it has been shown that factoriza- 
tion operators are directly related to the supersymmetric charges [U El [7] 
introduced by Witten [SJ. 

On the other hand, constants of motion for a given physical problem allow 
to simplify the corresponding Hamiltonian by reducing its degrees of freedom. 
Moreover, such conserved quantities are directly associated with symmetry 
groups and compact and non-compact Lie algebras [9]. Such algebras play 
a central role to study many properties of quantum systems because, among 
other things, they are the basis for selection rules that forbid the existence of 
certain states and processes. There is not a unified way to obtain the explicit 
form of the generators of compact and non-compact algebras, but several to 
find them, as it is shown in [TOT ITT] . However, in references [121 131 EEH [15] 
attempts to systematize the construction of the su(l, 1) Lie algebra genera- 
tors from factorization methods have been reported. Recently, in a series of 
papers it has been shown that the Schrodinger factorization operators can 
be used to construct the su(l, 1) algebra generators for relativistic and non- 
relativistic central potential Hamiltonians [T61 [TTJ, [18] - These results enhance 
the importance of factorization methods in solving quantum systems. 

The three dimensional relativistic Kepler- Coulomb potential is one of the 
solvable problems in physics. This is due to the conservation of the total 
angular momentum and the Dirac and Lippmann- Johnson operators [19], 
which reduces to the Runge-Lenz vector in the non-relativistic limit. The 
first two are due to the existence of spin, and the latter gives account of the 
degeneracy in the eigenvalues of the Dirac operator of the energy spectrum. 
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Symmetries and SUSY QM are intimately related. In fact, it has been shown 
that supersymmetry is generated by the Lippmann- Johnson operator [15]. In 
reference [16] we have studied the relativistic Kepler-Coulomb problem from 
an algebraic approach by using the Schrodinger factorization to construct the 
su(l, 1) algebra generators. This problem admits to be treated in other ways: 
analytical [201 ED E21 E31 Ei E5] and factorization methods [HI EH EE], shape- 
invariance [22] , SUSY QM for the first- [201 E7J and second-order differential 
equations [2H], two-variable realizations of the su(l, 1) Lie algebra [29J and 
using the Biedenharn- Temple operator |30j . 

In [31] Joseph studied the 1/r potential in D + 1 dimensions by means of 
the self-adjoint operators. Recently, the energy spectrum and the eigenfunc- 
tions of this problem were obtained by solving the confluent hypergeometric 
equation [32l [33] . Moreover, in [31] the Johnson-Lippmann operator for this 
potential has been constructed and used to generate the SUSY charges. 

The Dirac equation with Coulomb-type vector and scalar potentials in 
three dimensions has been solved by using SUSY QM [35] and by an analyt- 
ical approach [55] . In [27], this problem has been solved in D-dimensions by 
constructing the angular eigenf unctions from the group theory and the ra- 
dial equations solutions were expressed in terms of confluent hypergeometric 
functions. 

In the present Letter, we study the Dirac equation with Coulomb-type 
vector and scalar potentials in D + 1 dimensions from an su(l, 1) alge- 
braic approach. We obtain the uncoupled second-order differential equa- 
tions for bound states satisfied by the radial components and by applying 
the Schrodinger factorization we construct the corresponding su(l, 1) algebra 
generators. We use the theory of unitary representations to obtain the en- 
ergy spectrum and the action of the su(l, 1) algebra generators on the radial 
eigenstates. Also, we find the Schrodinger and SUSY ground states. By par- 
ticularizing our results to the cases where there exists either scalar or vector 
potential, or where both potentials are equal, we show that our treatment 
successfully reproduce the results obtained by analytical techniques. Finally, 
we give some concluding remarks. 
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2 The relativistic Dirac equation in D + 1 di- 
mensions 



The Dirac equation in general dimensions is 

H^= \ y^a aPa + /3(m + V s (r))+V v (r)J* = i—, (1) 

with h — c — 1, m is the mass of the particle, p a = —id a = — igf-, 1 < 
a < D, V s and V„ are the spherically symmetric scalar and vector potentials, 
respectively. Moreover, a a and (3 satisfy the relations a a oib + ctb®a = 2<5 a ftl, 
(%a(3 + 0cx a = and a 2 = (3 2 = 1 [37]. In Z) spatial dimensions the orbital 
angular momentum operators L a b, the spinor operators S a b and the total 
angular momentum operators J a b are defined as L a b = —Lb a = ^ x a^ — 
ix b -£^, S ab = -S ba = i^f 2 - and J ab = L ab + S ab , respectively. Thus, L 2 = 

Ef<> a6 , S 2 = Za< b S 2 ab, J 2 = Ea< b JL with 1 < a < b < D. Hence, 
for a spherically symmetric potential, the total angular momentum operator 

J a b and the spin-orbit operator Kd = —f3 ^ J 2 — L 2 — S 2 + J commute 
with the Dirac Hamiltonian, if. For a given total angular momentum j, the 
eigenvalues of Kp are = ± (j + (D — 2) /2), where the minus sign is for 
aligned spin j — £ + |, and the plus sign is for unaligned spin £ — \ [37J. 

We propose the wave function of the Dirac Hamiltonian H to be of the 
form 

»tt«)=r-^ fg^^)^, (2, 

V ZG ^ ( r ) F /m(^D) / 

being G^\r) the radial functions, Y^f^) and Y^(f2£)) the hyperspheri- 
cal harmonic functions coupled with the total angular momentum quantum 
number j and E the energy. Thus, the Dirac equation leads to the radial 
equation 




V s - V v + m + E 
V v - V s + m - E zb. 



r 




(3) 



In this work, we consider the Coulomb-type scalar and vector potentials 
V v = —-, V s = — — , with a and a' positive constants. In the following 
sections we study separately the cases a^a' and a = a'. 
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3 Case a ^ a' 

By introducing the new variable 

p = (aE + a'm) rj ' \j\a' 2 - a 2 \, (4) 
equation can be written as 



a t KD - --y—^p- [k d + i aE+a , m ) t KD , (5 j 



where 



A+FS = sgn(a - «')^R^ (*» ~ ^^S) ( 6 ) 
d 7 , V\ 



\a' 2 — a 21 



A ± = ±j-~ 1 + ^ ^, (7) 

dp p 7 



F,S\ _ / k d + 1 -{a-a')\ (gP 



(8) 



,1 if a > a' 

s#n(a - a ) = <( (9) 
-1 it a < a 



a 2 . 



and 7 = a/ k d + a' 2 

From equations §5$) and (JH]) we obtain the uncoupled second-order differ- 
ential equations 



d 2 7(7±1) \a' 2 -a 2 \ 2^\a' 2 -a 2 ' 



_ fVf^ 1 ; , I" ~ " I _ F ~ « I p (i,2) 



(ip 2 p 2 7 2 p 

1/2 2 ^ + am 



sgn (a — a') — \ n D — 7 



7 2 \ V aE + a 



where the superscripts (1, 2) correspond to (+, — ) in the centrifugal term. 
Thus, for bound states (m 2 > E 2 ), equation (jTDJ) is written as 

(V^ + *V - 2yJ\a*- a 2\p^ F« = -7 (7 + 1) F«, (11) 

(-p 2 ^ + *V - 2^ 2 -a 2 |p) *S = -7 (7 - 1) F£), (12) 
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with 



£ 2 = sgn(a — a') 



a'E + am 



- 1 



(13) 



aE + a'm / 

It must be noted that equation ( TTTT) is formally obtained from equation ( ITS]) 

(2) 

by performing the change 7 — > 7 + I. Hence, by defining -?/> 7 = Fk d , it implies 
that oc ^7+1- Thus, the solution for the Dirac equation in spinorial form 

can be expressed as •„ = ( F ™\ = (V 



3.1 The sw(l, 1) Lie algebra and the energy spectrum 

In order to construct the su(l, 1) algebra generators, we apply the Schrodinger 
factorization [21 IB] to the left-hand side of equation f[T2j) . Thus, we propose 



p4z + a P + b ) ( _ P:r: + C P + / ) ^7 = A'V'tj 



(14) 



<ip J \ dp 

where a, b, c, / and g are constants to be determined. 

Expanding this expressions and comparing it with equation ( TT2|) we obtain 

a = c = ±^, / = l + 6 = ^'"g"" 2 ' and # = 6(6 + 1) - 7(7 - 1). Therefore, 
the differential equation satisfied by the lower component of the spinor <3>fc D 
is factorized as 

2 



(JS^TI)^, 



± 



1 



7 



7> 



where 



dp £ 



(15) 



(16) 



are the Schrodinger operators. 

Also, we define the new pair of operators 

d 



n± = =Fp— + z P - n 3 , 

dp 



where 



1 ( d 2 _ 2 7(7-1) 

^ 1 = -(-p^+ep+ • 



2? 



dp 2 



p 



V>7 



v, 



a' 



-A 



(17) 



18) 
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has been obtained from equation (|T2"|) . 

A direct calculation shows that operators n_|_ and n 3 close the su(l, 1) 
Lie algebra 

[n±, n 3 ] = Tn±, [n + , n_] = -2n 3 . (19) 

Therefore, from equations (fTTj) and ffTB"]) we find that the corresponding 
quadratic Casimir operator II 2 = — n + n_ + (n 3 ) 2 — n 3 , satisfies the eigenvalue 
equation 

rr> 7 = 7(7-1)^. (20) 

As we mentioned above, by performing the change 7 — > 7 + 1 to 
we obtain the upper component of the spinor t/^+i- Thus, the su(l, 1) 
algebra generators for this component are 



]_(_ , ,2 , 7(7 + 1) 

2£ \ P dp- 



r± = =FP^ + fp-r 3 . (22) 

In order to determine the properties of the operators Il 3 and IT-|-, we intro- 
duce the inner product on the Hilbert space spanned by the radial eigenfunc- 
tions for the Dirac equation with scalar and vector potentials in D dimensions 



(0, C) = / 0* (p) C (P) P^dp. (23) 



Thus, the operator Il 3 is Hermitian with respect to this scalar product. 
Moreover, using equations ( ITT]) and ( 1231 . it can be proved that operators U± 
are hermitian conjugates, n_|- = n^. 

The theory of unitary irreducible representations of the su(l, 1) Lie alge- 
bra has been studied in several works [381 EH] and it is based on the equations 

T 2 |^)=/i(/i + l)l/W, (24) 
T 3 \fjiv) = v\nv), (25) 

T±\n u) = [(y^n){y±n± 1)] 1/2 \pu±l), (26) 



where T 2 is the quadratic Casimir operator, v = p, + q + 1, q = 0,1,2,... 
and \i > — 1. From equation fl26l) it can be noted that T + (T_) is the raising 
(lowering) operator for v. 
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Thus, from equations (I2"0"j) and 

= 7 - !, 

z/ 7 = n 7 + 7 



and ([TS]) and fl2H, we find that 



a" 



(27) 
(28) 



0,1,2, 



respectively, with n 7 

Since the operators IT-|- leave fixed the quantum number /z 7 , equation 
fT27|) ensures that the values of 7 remain fixed. This is because the change 
v 1 — > z/ 7 ± 1 induced by the operators n ± on the basis vectors |/i 7 z/ 7 ) implies 
n 7 — )• n 7 ± 1. Thus, by setting |/i 7 z/ 7 ) — )■ ?/> 7 7 and from equations (12"6T) . f[2"Tl) 
and fl28l). we find that 



C. 



(29) 



with C'j™ 7 ' 7 ^ = [(n 7 + 7 =F 7 ± 1) (n 7 + 7 ± 7)]^ 2 a real number. 

By using the equations (TT3"j) and (I2"8"j) we find that the energy spectrum 
for the lower component of the spinor $ KD , ip™~" , is 



£ 7 /m = [-a'a ± C^/C + a 2 - (a') 2 ] / (a 2 + ( 2 ) 



(30) 



where ( = 7 + n 7 . 

If we perform the change 7 — )■ 7 + 1 in equation ( |29|) . we find that the 
action of the ladder operators T± on the functions VvTi 1 is 



, rc 7+ i 



^,(n 7+ i,7+l) ,n 7+ i±l 



(31) 



KD 



r±^ 7+ i — ^± yvh 

meanwhile the energy spectrum for the upper component of the spinor $ 
is 

E 1+l /m = [-a'a ± CV(C') 2 + « 2 - K) 2 ] / (« 2 + (0") , (32) 

where £' = 7 + ^7+1 + 1- However, since ^ 7 +V and ?/> 7 7 are the components 
of the same spinor, they should have the same energy, E 7 = E 1+ i. Thus, 
we obtain n = n 7 = n 7+1 + 1, where n = 0, 1, 2, 3, ... is the radial quantum 
number. 

Therefore, the energy spectrum for the Dirac Hamiltonian with scalar 
and vector potentials in D dimensions can be written as 



E 



m 



-a'a 



± 



a 2 + (7 + n) 2 ^ \a 2 + (7 + n) 2 / \ « 2 + (7 + n) 



a 



a 



(7 + n\ 



(33) 
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which is the same reported in [37J , obtained from an analytical point of view. 
Moreover, this expression is reduced to that for the three dimensional case 
[36]. 

In this way we have shown that $™ is given by 

^ = (it) = (If) ■ (34) 

The relation between the components of the spinor above via the radial 
quantum number n, deduced from the theory of unitary representation can 
be observed from an analytical approach, as it is shown below. 

3.2 The Schrodinger and SUSY QM ground states 

From equation ( 12 9 p . for n = we find that only the state ip® SCH = r^e""^' ' 2- " 2 '^ 7 
is normalizable with respect to the inner product defined in expression ( 1231 
and satisfies the differential equation IT_?/>° = 0. Moreover, for n = and 
from the definition of the radial quantum number, we find that n 7+1 = — 1. 
For this value, C_ results in a complex number. From the theory of uni- 
tary representation and equation ( 13T1) . the function VC+i is non-normalizable 
[39] . Therefore, this function is not a physically acceptable solution and the 
spinor corresponding to n = is 



$ 7 SCH - ( ^y e -V|a' 2 -" 2 |p/7 I ' ^ 



which we denote as the Schrodinger ground state. 

In fact, the operators A ± in equations (jSJ) and (EJ) are the SUSY operators, 
from which the partner Hamiltonians are given by H + = A~A + and H = 
A + A~. The ground state for H_ satisfies the condition ^4 _ ^susy = 0' tlmt 

leads to the square-integrable eigenfunction SUSY = cp 7 exp (— p^\a' 2 — a 2 \/^fj . 

On the other hand, the solution for the equation A + ip® +1 SUSY = is ?p® +1 SU sy = 

cp" 7 exp ^— pa/|«' 2 — a 2 |/7^. Since this function is not a square-integrable 

then, it is not a physically acceptable solution and it must be taken as the 
zero function. Hence, the SUSY ground-state is given by 

$ ?SUSY = L 7e -V|a' 2 -a 2 |p/ 7 J ' ( 36 ) 
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The above results leads us to the conclusion that the Schrodinger and SUSY 
ground states for the relativistic Dirac equation with Coulomb-type scalar 
and vector potentials are the same. 

The explicit form of corresponding to higher energy levels can be 
obtained from an analytical approach. In order to solve the differential equa- 
tions (TT21) we propose F^j = p 7 e~^ p /(p), where f(p) must satisfy 



v ,J 2 -a 2 \ 

y-r^ + K 2 ! -y)-r + ? 7 



d 2 d vT"" ~~ n ~ 



dy 2 dy £ 



/(j//2O = 0, (37) 



with y — 2£p. The solution for this equation is the confluent hypergeometric 
function 1 F 1 (-n,2j;2£p). Thus, we get F r ?J D = r] 2 p 7 e"^ x Fi (-ra, 27; 2£p). 
In a similar way, the solution for equation ffTTl) is FnJ D = r]ip 1+l e~^ p 
±Fi (— n + 1, 27 + 2; 2£p). The relation between these solutions is in agree- 
ment with that obtained from the theory of unitary representations, equation 
( |34|) . Moreover, it is known that i-Fi(0, b; z) = 1, whereas iF±(a, b; z) diverges 



for a > 0. Therefore, for n = 0, Fq^ = rj 1 p' y e~^ p , while the function Fq^ 
is not square-integrable, and it must be taken as the zero function. This 
result is in accordance with equation ( 135]) . which has been obtained from an 
algebraic approach. 

From the definition of n, equations fT29|) and fl3T|) imply 



n ± ^ocVf\ r^^oc^r 1 . (38) 



It is worth to notice that equation f[38|) is valid for any value of the radial 
quantum number except for n = 0. This is because the upper component of 
the spinor for n = 1 can not be obtained from the action of the operator T + 
on the upper component of the Schrodinger ground state, equation (1351) . 

Equations (fT8|) and (|38|) allow us to show that the action of the Schrodinger 
operators on the states ip™ is ££±i\?l oc , 0~ ±1 . A similar result can be ob- 
tained from the operators r±. Therefore, we conclude that the action of the 
Schrodinger operators Jzf± on the components of the spinor is to change 
only the radial quantum number n leaving fixed the Dirac quantum number 

«D = «I>(t)- 

Some particular cases of our general results obtained above are shown in 
table 1. 

From an analytical point of view, case I was studied in [32J, meanwhile 



cases III, IV and V were treated in |36j. On the other hand, case III was 
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Cases 


Spatial 
Dimensions 


Potential 
Parameters 


7 


e 2 


/ 


D 


a ^ 0, a ' = 




(f)'-l 


II 


D 


a' ^ 0, a = 






III 


3 




v/«3 + a' 2 — a 2 




IV 


3 


a ^ 0, a ' = 


a/k| - a 2 


(tr-1 


V 


3 


a' ± 0,q = 




!-(-)" 

V m / 



Table 1: It is shown the expressions for 7 and £ 2 for particular values of a and a 1 
in L> or three spatial dimensions. 



studied by SUSY QM in [35] and case IV was treated algebraically in [16J. 
It must be pointed out that our present results, which were obtained from 
an su(l, 1) algebraic approach, are in full agreement with those reported in 
the last references. 



4 Case a = ol 

Because of the change of variable given in , a = a' is not an special case of 
the results of the previous section. For the present case, from the Dirac radial 
equation (j3J), we obtain the uncoupled second-order differential equations 

(-P 2 |^ + eV - %*p\ *S = -Kd {kd + 1) F«, (39) 
(-p 2 ^ + £V - 2ap) F£> = -k d {k d - 1) F#, (40) 
where £ 2 = ^r§, p = {m + E)r and 




(41) 



Since equation ( 1391) is formally obtained from equation ( 1401) by performing 
the change kd — )■ + 1 then, by defining ijj KD = F« D , we conclude that the 
solution to the Dirac equation in spinorial form is 



•» s 1 1) = Ct +1 ) ■ « 42) 
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Proceeding in a similar way as in the case a 7^ a', we obtain that the 
su(l, 1) algebra generators for the upper and lower components of $ Kc are 



X3 = X3[K D j 



2£ V P dp 
d 



and 



X± = X±(«z>) = TP^ +£p~X3 



A3 = X3<^D + 1), A± = X±(«£» + 1), 



(43) 
(44) 

(45) 



respectively, with X3^k d = WOV^d and A 3 ^ KD +i = (a/OVVo+i- 

From the theory of unitary representations, equations ff24|) . (125]) and (T56 
we obtain that the corresponding energy spectrum is 



E = m 



2a 2 



a 2 + (n + ko 



(46) 



Also, we find that the Schrodinger and SUSY ground states are the same 
and are given by 



^SCH,SUSY 





p K D g-ap/ftu 
l K £)| e - a P/l K D 





for kd > 0, 
for kd < 0. 



(47) 



These results are in accordance to those reported in [35] for three spatial 
dimensions. 



5 Concluding remarks 

We studied the radial Dirac equations with Coulomb-type scalar and vector 
potentials in D + 1 dimensions from an su(l, 1) algebraic approach. With 
the Schrodinger factorization we were able to construct two sets of su(l, 1) 
algebra generators. We applied the theory of unitary representations for this 
algebra to find the general form of the spinor wave function, the energy spec- 
trum and the SUSY ground state. We proved that Schrodinger and SUSY 
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ground states are the same. We showed that the action of the Schrodinger op- 
erators on the radial eigenstates is to change only the radial quantum number 
n leaving fixed the Dirac quantum number Kd- The noncompactness of the 
su(l, 1) algebra reflects that, for a fixed number kd, the quantum number n 
is bounded from below and unbounded from above. We tested our results by 
matching different cases where there exists either scalar or vector potentials 
in D or three spatial dimensions [TBI E21 ESI EH EI]- To our knowledge this 
is the first time where the Dirac equation in D dimensions with both scalar 
and vector potential has been treated by an su(l,l) algebraic approach. 

Finally, our technique can be successfully applied to solve other relativis- 
tic problems like the Coulomb field with position-dependent mass [40] or the 
Dirac- Morse [4T], which are works in progress. 

Acknowledgments 

This work was partially supported by SNI-Mexico, COFAA-IPN, EDI-IPN, 
SIP-IPN project number 20110127 , and ADI-UACM project number 7DA2023001. 

References 

[1] Dirac P. A. M., The Principles of Quantum Mechanics, Clarendon Press, 
Oxford, 1935. 

[2] Schrodinger E., Proc. R. Ir. Acad. A 46(1940)9. 

[3] Infeld L. and Hull T. E., Rev. Mod. Phys. 23(1951)21. 

[4] Cooper F., Khare A. and Sukhatme U., Phys. Rep. 251(1995)267. 

[5] Dutt R., Khare A. and Sukatme U. P., Am. J. Phys. 56(1987)163. 

[6] Bagchi B. K., Supersymmetry in Quantum and Classical Mechanics, 
Chapman & Hall/CRC, USA, 2001. 

[7] Lahiri A., Roy P. K. and Bagchi B., Int. J. Mod. Phys. A 5(1990)1383. 

[8] Witten E., Nucl. Phys. B185(1981)513. 

[9] Wybourne B. C, Classical Groups for Physicists, Wiley Interscience, 
New York, 1974. 



13 



[10] Bohm A., Ne'eman Y. and Barut A. O., Dynamical Groups and Spec- 
trum Generating Algebras, vols. 1 and 2, World Scientific, Singapore, 
1988. 

[11] Englefield M. J., Group Theory and the Coulomb Problem, Wiley In- 
terscience, USA, 1972. 

[12] Miller W., Lie Theory and Special Functions, Academic Press, New- 
York, 1968. 

[13] Levai G., J. Phys. A: Math. Gen. 27(1994)3809. 

[14] Gangopaghyaya A., Mallow J. V. and Sukhatme U. P., Phys. Rev. A 
58(1998)4287. 

[15] Martinez D., Granados V. D. and Mota R. D., Phys. Lett. A 
350(2006)31. 

[16] Salazar-Ramirez M., Martinez D., Mota R. D. and Granados V. D., J. 
Phys. A: Math. Theor. 43(2010)445203. 

[17] Martinez D., Flores-Urbina J. C, Mota R. D. and Granados V. D., J. 
Phys. A: Math. Theor. 43(2010)135201. 

[18] Martinez D. and Mota R. D., Ann. Phys. 323(2008)1024. 

[19] Dahl J. P. and Jorgensen T., Int. J. Quant. Chem., 53(1995)161. 

[20] Thaller B., The Dirac Equation, Springer- Verlag, Berlin, 1992. 

[21] Thaller B., Advanced Visual Quantum Mechanics, Springer- Verlag, 
Berlin, 2005. 

[22] Swainson R. A. and Drake G. W. F., J. Phys. A: Math. Gen. 24(1991)79. 

[23] Swainson R. A. and Drake G. W. F., J. Phys. A: Math. Gen. 24(1991)95. 

[24] Su J., Phys. Rev. A 32(1985)3251. 

[25] Wong M. K. F., Phys. Rev. A 34(1986)1559. 

[26] de Lima Rodrigues R., Phys. Lett. A 326(2004)42. 



14 



[27] Sukumar C. V., J. Phys. A: Math. Gen. 18(1985)L697. 

[28] Jarvis P. D. and Stedman G. E., J. Phys. A: Math. Gen. 19(1986)1373. 

[29] Martmez-y- Romero R. P., Nunez- Yepez H. N. and Salas-Brito A. L., 
Phys. Lett. A 339(2005)259. 

[30] Horvathy P. A., Rev. Math. Phys. 18(2006)311. 

[31] Joseph A., Rev. Mod. Phys. 39, (1967)829. 

[32] Jiang Y., J. Phys. A: Math. Gen. 38, 1157(2005). 

[33] Dong S. H. and Lozada-Cassou M., Int. J. Mod. Phys. E 5(2004)917. 

[34] Katsura H. and Aoki H., J. Math. Phys. 47(2006)032301. 

[35] Guo-Xing J. and Zhong-Zhou R., Commun. Theor. Phys. 49(2008)319. 

[36] Greiner W., B. Miiller and J. Rafelski, Quantum Electrodynamics of 
Strong Fields, Springer- Verlag, Berlin, 1985. 

[37] Dong S. H., Sun G. H. and Popov D., J. Math. Phys. 44(2003)4467. 

[38] Adams B. G., Cizek J. and Paldus J., Adv. Quant. Chem. Vol. 19, 
1(1987). Reprinted in A. Bohm, Y. Ne'eman, A. O. Barut, Dynami- 
cal Groups and Spectrum Generating Algebras, vols. 1 and 2, World 
Scientific, Singapore, 1988. 

[39] Adams B. G., Algebraic Approach to Simple Quantum Systems, 
Springer- Verlag, Berlin, 1994. 

[40] Alhaidari A.D., Phys. Lett. A. 322(2004)72. 

[41] Alhaidari A.D., Phys. Rev. Lett. 87(2001)210405-1. 



15 



